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Abstract 

We investigate the effects of a Stern-Gerlach-type addition to the Lorentz force on electrons in a laser 
wakefield accelerator. The Stern-Gerlach-type terms are found to generate a family of trajectories describing 
electrons that ‘surf’ along the plasma density wave driven by a laser pulse. Such trajectories could lead to an 
increase in the size of an electron bunch, which may have implications for attempts to exploit such bunches 
in future free electron lasers. 


1 Introduction 

The modelling of charged particles in electromagnetic fields has been an area of great interest to accelerator 
physicists for decades, using a number of simplifications for computational ease. Increasingly large fields such 
as those that will be produced at ELI [T] and HiPER [5] now mean that the validity of these simplifications 
must be considered more closely. 

The classical radiation reaction force has long been considered as an averaged effect for sufficiently low 
electromagnetic fields in accelerator physics and more recently the Landau-Lifshitz equation has been used to 
model the motion of such particles in detail [aiiiis]. The Landau-Lifshitz equation is however only valid for 
sufficiently weak and slowly varying background fields [5] . The increasingly large electromagnetic fields expected 
to be deployed at experiments such as ELI and HiPER may be strong enough to render the Landau-Lifshitz 
equation invalid and this has led to recent interest in alternative models [7], computational simplifications [5] 
and extensive consideration of quantum effects (for a recent review, see i)- Thus, calculation of the effective 
trajectories from first principles proves problematic, due to the various candidate effective theories. Recently 
[8], various radiation reaction models have been been included in PIC codes via the Osiris framework, which is 
commonly used in the laser wakefield context. 

While the effects of the increasingly large electromagnetic field strengths have been studied, the effects of 
large field gradients, created in laser wakefield accelerators, have received very little attention. In contexts 
where the field gradient is considered to be non-trivial the radiation reaction contributions are included via the 
Landau-Lifshitz equation, though the effects of spin are often neglected despite being of comparable size [9] . The 
Stern-Gerlach force, caused by the coupling of the spin of a charged particle with the background electromagnetic 
field, was first observed in the splitting of a beam of silver atoms by an inhomogeneous magnetic field. In high 
field-gradient systems, such as those created in ultra-intense laser-driven plasma density waves, we suggest that 
the Stern-Gerlach force is not only non-trivial but in some circumstances may be more important than the 
radiation reaction force. These Stern-Gerlach effects may lead to behaviour that does not appear to have been 
adequately addressed in the context of laser wakefield accelerators: the purpose of this article is to illustrate 
the significance of the Stern-Gerlach forces in a simple model of such an accelerator. 

Modelling quantum mechanical electrons as covariant classical spinning particles has been well-studied. 
There have been various approaches from the work of Erenkel m and Thomas m in the 1920s through the 
work of Nakano [12], Tulczyjew [13], Dixon [HIIlllIl], Gorben pTl [T8] , Suttorp and de Groot [TU |20] and Ellis 
[2T] in the 1950-70s. The equations of motion and the validity of the auxiliary condition on the spin and the 
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momentum necessary to close the system have still received discussion in recent years [22] . The approaches used 
to derive these equations are varied, and for completeness we include a new method using de Rham currents 
and distributional methods in Appendix [S] 

The systems in which the Stern-Gerlach force is most prominent are those with a high electromagnetic field 
gradient. Chapter 2 considers the implications of the coupling between the spin of a classical electron and the 
rapidly varying electromagnetic field produced by a laser-driven plasma wave. Sufficiently short, high-intensity 
laser pulses can form longitudinal waves within the electron density of a plasma. These density waves propagate 
with speed comparable to the group speed of the laser pulse. Not all plasma electrons form this wave, however; 
some of the electrons are caught up in the wave and accelerated by its high helds. The wave eventually collapses 
as these electrons damp the wave (the wave ‘breaks’). The extremely high electric field gradient of a plasma 
wave near wavebreaking provides an excellent theoretical testing ground for the effects of Stern-Gerlach-type 
contributions to the trajectory of a test electron. 

In what follows, the equations of motion for a spinning electron in such a density wave are found to have 
a particular solution which does not exist for a particle without spin - trajectories corresponding to electrons 
‘surfing’ orthogonal to the wave vector in the frame of the wave. The perturbations around a ‘surfing’ trajectory 
are found to be linearly unstable for the vast majority of the parameter space. Since the family of new trajectories 
found in Section 2 correspond to electrons travelling orthogonal to the motion of the plasma electrons and are 
unstable, the electrons following such trajectories could cause undesirable properties for effective bunching of 
electrons in laser wakefield accelerators. 

These ‘surfing’ trajectories exist only for a particle with non-zero spin in a background held with non¬ 
zero gradient. Furthermore, the electrons are non-accelerating and therefore the radiation reaction forces are 
expected to be negligible. Clearly, the spin-held coupling is much more signihcant than radiation reaction in 
the present context. 


We use Heaviside-Lorentz units with the speed of light c = 1 (except at the end of Section 2.3 for the sake 
of clarity) and we assume that the effects of spacetime curvature are negligible so that the spacetime metric is 
simply the Minkowski metric gat = diag{—1,1,1,1}. Lower case Latin indices run over 0, 1, 2, 3. 


2 Effects of Stern-Gerlach-type forces on a classical charged particle 

2.1 Preliminaries 


The equations of motion that govern a classical particle with worldline C : r i—)■ a;“ = C“(r), charge q, momentum 
P“ and spin Sab in a background electromagnetic held described by the tensor components Fab are 


^ , F-^i:bcP^ \ 

dr \ Cdp^ ) 


-qF^^Cb - 


±s^b ^ fpb ^ ^ fp. p 

dr \ C^Pa J V C-Pe J 






( 1 ) 

( 2 ) 


where r is the proper time of the particle, (7“ = is the 4-velocity of the particle and Y,ah is the 

electromagnetic dipole tensor (see |19j . or for a new derivation of these equations, see Appendix A). However, 
this is not a complete system; an additional condition is required. There are a number of possible conditions, 
though two of the most commonly used are the Frenkel condition m 

C^Sab = 0 (3) 


and the Nakano-Tulczyjew [HIS] condition 


P'^Sab = 0 . 


(4) 


The Frenkel condition, whilst being simple and intuitive, is considered by some to be unphysical since it yields 
helical solutions in held-free systems (sometimes called Zittebewegung) [T71 [TH] , though others argue against this 
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unphysicality [^. We initially adopt the Nakano-Tulczyjew condition Q, which has already been abundantly 
studied [H EH m ng though we subsequently show that in fact the conditions <§ and @ are equivalent to 
first order in Sab- 

A particle with spin has a magnetic dipole moment related to the spin by the gyromagnetic ratio , where 
Mq is the particle’s rest mass and g is the g-factor of the particle. Furthermore, a particle with a zero electric 
dipole moment and a non-zero magnetic dipole moment is characterised by Eab = ^^^Sab- Thus the equations 
of motion, together with the condition Q, are 


dpa ^ _ ^s'^^d^^Fba, 

dr 4Mo 

gab _ _ pb i ^bpa i 9Q pbcg _ 99 pacg b 

dr ~ 2Mo 2 Mo 


(5) 

( 6 ) 


where the 4-momentum satisfies the condition 


Pb 


-P'^CaCb - 


{Sab + SadC'^Cb 


99 


9F’^^Ca + i^d^FadS^^^ + 

4iIVIr) 


99 

2Mn 




(7) 


found by differentiation of the Nakano-Tulczyjew condition Q with respect to t. Note that the first term on 
the right-hand side of ([^ is the standard Lorentz force on a charged particle and the second term, the coupling 
of the spin and the gradient of the electromagnetic field, is a Stern-Gerlach-type contribution. 

A classical electron has g-factor equal to 2, charge q — qe = —e (where e is the elementary charge) and rest 
mass Mq = me giving the system of equations 


= -qeF-^Cb - ^S^^d^Fba, 
dr Zrrie 

^ gf^^ -j- P^ F pbcg a 9e pacg b 

dr me me ’ 

P^Sab = 0 , 

pa ^ _pb(j^^ca _ 


( 5 ,“ + SbdC'^C°) 


P^Ce 


qeF^^^Cc + —d'^FedS^'^ + -^P"P 

rrif. TTif. 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Integration of is far from straightforward due to constraints ( [l 0 | ) , 0 and it is useful to reduce ([H)-([Tl|) 

to a model that captures the essential physics that we wish to explore. In order to simplify the system, we 
choose to linearise (|H)-(fll|) in the spirQ S'ah- Firstly, note that linearising the momentum condition (II) results 
in the straightforward expression 


for the 4-momentum and hence the Frenkel condition ([^ and the Nakano-Tulczyjew condition (|^ are equivalent 
to first order in 5°“^. The linearised system of equations is 




9e 

2 to ? ‘ 


gbcpb^ 1 ^F'^^Cb - 


9e 


<ab/ 


^ab _ S ^ _ F^^S 

rln- rKt ' ^ 


o = 

dr m, 

& Sab = 0 


me 

b\ 


9e 

2m‘i 


s'^^ {n^y^'ddFb^ 


(13) 

(14) 

(15) 


along with (12), where (11^)“° = 5 “° -|- C°'C°. Notably, the effects of the third term (the Stern-Gerlach-type 


term) on the right-hand side of (13) are most apparent in a system with a high field gradient. Although we 


^For a system of equations linearised in Fab-) see m- 
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can demand that (151 is satisfied at a particular instant in proper time r, it is only satisfied to second order in 
5“^ at other times. Henceforth we adopt (13), (14) as the system of equations for a classical electron with spin 
subject to C°'Sab\T=Q = 0 in order to readily demonstrate the effects of the Stern-Gerlach term (the final term 
on the right-hand side) of (13). 


2.2 Effects of the Stern-Gerlach force on the motion of an electron in a plasma 
wave 

For clarity we use the simplified notation t = (7°, x = y = C^, z = for the components of the worldline 
C. Consider a system with the electromagnetic field 2-form F associated with an electrostatic wave 


Fah = 


E{^) for a = 0, 6 = 3 
—E{^) for a = 3, 6 = 0 
0 otherwise, 


where ^ = z — vt is the phase of the wave. The evolution equations (13) are simply 


i = f 1 + -^SosE^ —Ez + -^ {v + i{z- vi)) E'Sqs, 

\ J rrig Wg 

X = -^x (i - vi) E'Sos, 

y = {z - vi) E'So3, 

z = (l + -^SosE^ —Ei + -^ {l + z {z - vi)) E'Sqs, 

\ rrig J rUe 


(16) 


(17) 

(18) 

(19) 

( 20 ) 


since = 


C and Fo 3 = —F^^. Note that dots denote derivatives with respect to the proper time r and 


primes denote derivatives with respect to the phase Similarly from (14) the spin evolution equations are 


bi 

= —ES 33 , 

^13 = 

— ESoi, 

( 21 ) 


me 


rUe 

^2 

= —ES 23 , 

S 23 = 

—ES 02 , 

( 22 ) 


me 


me 


= 0, 

S 12 = 

0 . 

(23) 


Notably, the only component of the spin in @-([ 2 ^ i.e. the only component that affects the trajectory of the 
particle is S' 03 , which according to (23) is constant. We hence neglect the remaining spin equations of motion 
when solving for the worldline of the electron. Writing the remaining equations of the system, ([^-(l^, in 
the coordinate system { 7 C; 2 ;, y, 7 ^}, adapted to an observer travelling with the plasma wave at speed v, where 
7 = , ^ is the Lorentz factor of the wave and C = —t + vz, we find 


C 

I 

X 

y 


^So3E'a- 

mi 

^So3E'U- 

mi 

^So3E'xt 

mi 

^So3E'yi 

mi 


(1 + ||S.3£) 


qe 

me 

qe 

me 


Ei 


EC + 


qe 

mi^"^ 


S 33 E' 


(24) 

(25) 

(26) 
(27) 
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A particular solution to (24)-(27) includes a constant value for phase ^ and has the form 


Csoi(t) — 


*S'o3 Eq 


(l + Ec 

Xso\{t) = XqT + Xq, 

ysoi(T) = yoT + yo, 

^sol = Cc, 


T + Co, 


(28) 

(29) 

(30) 

(31) 


and where Eq denotes the value of the electric field at ^ = ^c- Here Xq, yo, Co are arbitrary constants and xq, 
yo are arbitrary constants up to fulfilment of the normalisation condition gabC°‘C^ = — 1 on the worldline of the 
electron, i.e. 


(i+^503Ac)^ 


So?, E'r, 


TOe7' 


03 

2 Ac 


yl = -1- 


(32) 


Note that condition (32) places restrictions on the system parameters, for instance Sq?, Eq, Eq ^ 0. Con¬ 
sequently, the solution family (28)-(31| does not exist for a spinless particle, nor a constant electromagnetic 


ft to 


background. Also note that since we wish to consider a system with a large field gradient, we consider 
be of order (S'os)”^. 

A sufficiently short and intense laser pulse propagating through a plasma may create a travelling longitudinal 
plasma wave whose velocity is approximately the same as the laser pulse’s group velocity. The electric field 


produced by such a plasma wave provides an excellent example of an electric field of the form (16): 


E = 


QeT 


(33) 


where v must satisfy 


^2^ ~ ^eZriion (u-v/l/2 - ^2 

2qij^ V 



(34) 


in order to satisfy the Maxwell equations and the Lorentz force equation (see Appendix for details). Here 
Z is the degree of ionisation, riion is the proper number density of the (background) ions. In laser wakefield 
acceleration, the “target” is the dephasing point, where accelerated electrons begin to overtake the plasma wave. 
At this point the field gradient is much larger than the electric field, hence ^ can be said to be large (of order 

(^03)-^). 

The family of trajectories given by (28)-(31) with the electric field (|33|) is illustrated in Figurej^ Despite the 


propagation of the plasma electrons in the ^ direction, the electrons described by the solution family (|28|)-(|31|) 
travel transversely, along lines of constant E (they ‘surf’ along the wave). 


2.3 Stability of the ‘surfing’ solutions in a plasma wave 


It is natural now to consider the linear stability of the family of ‘surfing’ solutions described by (28)-(31) for 
the plasma wave electric field (33). In order to investigate this, consider the following: 


C(t) = Csoi(r) -h£AC(T), 

(35) 

x{t) = Xsoi(t) -I- eAxir), 

(36) 

yix) =ysoi(T) +sAy{T), 

(37) 

((x) = ^so1+£A^{t), 

(38) 
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Figure 1: Illustration of several example trajectories Ci, C 2 , C 3 , given by different choices of ^c- Whilst the 


plasma electrons travel along test electrons described by (28)-(311 travel transversely to the wave’s velocity, 
‘surfing’ along the wave. 


where e is a small constant and the A terms correspond to perturbations. Substituting (38) into v and taking 
Taylor series in e gives: 


+ eA^) = vc + + 0{e^), 


(39) 






and so on for the higher derivatives. Substituting (39), its derivatives and 


where vc = v{^c), v'q = 

(35)-(38) into the equations of motion (|24|)-(|27[) gives, to first order in the perturbations. 


AC = 


-(1-5^^)^ + 
T 


52 


i^'6) 




1-Siy'r 


A^ = A AC, 


Ax = [Sxov'c] AC = ^ 2 AC, 

= [Syoi^c] AC = y^sAC, 

AC=^^i c„/^2(^c)V 


(1 - Su'c) 


S 




siiy'6) 


//^2 




AC = ^4AC, 


(40) 

(41) 

(42) 

(43) 


where S = — and the constants An depend on the spin, electric field and the plasma wave speed. A 
solution to (|40)-([43|) is 


AC = ^ +C2e-^") , 

Ax = ^ +C 2 e-^") , 

Ay = ^ +C 2 e--^") , 

AC = Cie^^-k C2e-^^ 


(44) 

(45) 

(46) 

(47) 


where Ci, C 2 are integration constants. The stability of the system hence depends solely on the sign of the 
quantity A 4 , defined in (43). Written as a Taylor series in S, A 4 can be expressed as 


.44 





+ 0{S^). 


(48) 
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Assuming that the zeroth and first order terms are dominant we neglect the higher order term^ Consequently, 
the exponential terms in the perturbations (44)-(47) become 


hcl. 


to first order in S, where 


gVMj. ~ e r 




(49) 


(50) 


Hence the ^ perturbation (38) to first order in S is 

beEcI . 


= Cl (1 + A^c^r) e + C 2 (1 — AqiSt) e’ 


beEcI . 


(51) 


that is the perturbation is unstable (to first order in 5) as the first exponential will diverge as r increases, 
unless Cl = 0. Since the other three perturbations are closely linked to A^, the complete perturbation is also 
divergent (unless the integration constant Ci = 0 ). 


2.4 The range of validity of (49) 

In order to confirm that the zeroth and first order terms of S in A 4 are dominant, consider the ratio R of the 
sum of the zeroth and first order terms to the full expression for A 4 (43). The ratio R can be written as 


R = 


(1 - A) 


- 2j^A')^ + A f - {AY'^ 


{y'Y + ~ 'y^A'^A 


(52) 


where A = 

2 3 

Returning to SI units for clarity, introducing the Schwinger limit and the maximum electric field 

(the wave breaking limit) for a cold plasma, i?max = allows the convenient re-parameterisation 

of the system in terms of the free parameters {f, A, k} where 


A = 

k = 


E So3 

Emax \ / <5'o3 \ 


(53) 

(54) 


The Schwinger limit characterises the electric field strength at which quantum vacuum effects are expected to 
be significant [5j . The present analysis does not take such phenomena into account and therefore the conditions 
E < Es, Ejnax < Es are required. Thus, in addition to 0 < r; < 1, it follows that —1 < < 1 and 0 < fc < 1. 

We also have the restriction 


k 



(55) 


to ensure that E < The combination of these conditions also guarantees that 1 / > 7 , ensuring that the 

square root in (34) remains real. The maximum amplitude plasma equation (341 and its derivatives relate A 


^The range of validity of this assumption is ascertained in the subsequent section. 
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Figure 2: The relative size of terms R in the parameter space {v’, k) for four values of plasma wave speed 
V. For each speed there are two plots, the upper using (56) with the positive sign and the lower using the 
negative sign. The regions in which the trajectories of ‘surfing’ particles (|2^-(31 ) are linearly unstable (i.e. the 
assumptions made in Section 2.3 are valid) are the regions where \R\ 1. 


and its derivatives: 


D± = -7^ ( (7 - 1)^^ - 7 I ± 7^1^ 


\ 


(7- 1) 




- 1 , 


1 k 




27-1 


- 1 


4-72 


(56) 

(57) 

(58) 


Figure illustrates the size of (52) across the parameter space {v', k) - it is clear that |i?| 1 for the 

majority of the parameter space. Several things should be made clear, however. Firstly, the black region in 
each plot is excluded by the condition (55) - these dark regions correspond to electric fields E > Flmax- The 


central line present in some plots indicates that the electric field must be non-zero, as is already stipulated by 


the normalisation condition (32). Secondly, the result of the numerical analysis is unreliable in certain regions 


of the plots, for example near the edge of the parabola or along the line v' = 0. 

For values of the electric held E of the plasma wave (and the maximum electric held of the cold plasma 
7?max) that are much less than the Schwinger limit E^, both \ v'\ and k are much less than 1 - hence the crucial 
regions of the parameter space to consider are relatively close to the origin. Figure [^illustrates the size of (52) 
across a smaller range of electric helds. Here in almost every plot the relative size of terms (52) is typically much 


less than 1, and hence the assumption (49) is generally valid. Hence it is safe to conclude that the trajectories 


described by (28l-(31) corresponding to the ‘surhng’ particles are linearly unstable across the vast majority of 


the parameter space for E Es- 


















































Figure 3: The relative size of terms R in the parameter space {v', k) for four values of plasma wave speed 
V. For each speed there are two plots, the upper using (56) with the positive sign and the lower using the 
negative sign. The regions in which the trajectories (28)-(31) of ‘surfing’ particles are linearly unstable (i.e. the 
assumptions made in Section 2.3 are valid) are the regions where \R\ <C 1. 


3 Conclusion 


A new family of trajectories for a classical charged particle with spin in an electrostatic plasma wave has been 
presented - notably these trajectories do not exist for a non-spinning particle, nor for a non-varying electric 
field. These trajectories represent particles moving transverse to the wave propagation, ‘surfing’ along the wave. 
The linear stability of these trajectories depends on the values of the plasma wave speed, the electric field and 
the spin component S 03 . As shown in Section |2.3[ these trajectories are linearly unstable for the majority of 
the parameter space, though there are some small regions in the parameter space where this may not be the 
case, where the assumptions of (49) are invalid. For lower electric field, as in Figure]^ these regions are even 
less prominent. 

The existence of transverse trajectories has adverse consequences for the size of electron bunches in laser 
wakeheld accelerators; electrons may slip into and out of these transverse trajectories once they catch up with 
the wave, spreading out into a disc oriented with normal parallel to the wave propagation direction. The linear 
instability of these solutions, however, implies that any electron that enters a transverse trajectory would likely 
leave it soon afterwards. 

It is important to note, however, that this is an instance in which the spin of a particle affects its trajectory 
more than radiation reaction effects - indeed in this instance, since the particles in the new trajectories are 
travelling at constant speed, the radiation reaction effects are negligible. Hence we recommend that researchers 
aiming to model such systems consider including spin effects in their PIC codes. 


4 Acknowledgements 

This work was undertaken as part of the ALPHA-X consortium funded under EPSRC grant EP/J0I817I/1 and 
with support from the Cockcroft Institute of Accelerator Science and Technology (STFC grant ST/G008248/1). 
S.P. Flood was funded by an EPSRC studentship. The authors thank Adam Noble for helpful comments on 


9 












































the manuscript. 


A Appendix: Derivation of the equations of motion 

This Appendix makes use of Euclidean 3-vectors, 4,-vectors on Minkowski spacetime and de Rham currents 
(Schwarz distributions on differential forms). For clarity, a different notation will be used for each of those 
objects. Use of an arrow V indicates a Euclidean 3-vector, V is a 4-vector with the appropriate metric dual V, 
and a subscript D such as Vo indicates a de Rham current. 

This Appendix shows a new derivation of the relativistic Stern-Gerlach force and TBMT equations for a 
charged particle, using exterior calculus and Schwarz distribution^ The aim of this section is to reduce the 
local balance laws [2l] 

= ix‘^FA + ix-^FA (59) 

i (da;“ AT'' - dx^ A r“) , (60) 

for a charged continuum, with stress 3-forms and spin 3-forms cr“^, to a particle model using de Rham 
currents. The vector basis {A“} is a Killing frame, and the electromagnetic 2-form F and current 3-forms 
jbound gg^|;jgfy Maxwoll oquatious 


dF = 0, (61) 

d*F=/"®®-hjb™''b^ ( 02 ) 

Given a charged fluid, the magnetisation and polarisation vectors are given by 

V{r,t) = n{r,t)fle{f,t), (63) 

M(f,t) = n{r,t)fL^{r,t), (64) 


where n is the particle number density and fie and /Xm are the electric and magnetic dipole moments respectively. 

In order to efficiently move from the continuum model to a single-particle model, de Rham currents [25] are 
introduced. Firstly, in order to establish the notation in a simple setting, it is assumed that the fluid is at rest 
and hence described by a 4-vector field given by V = 9t- Then the distributional current associated with the 
worldline of a particle is introduced via 


/n * 1 ^ / fdt, 


/M 


(65) 


analogous identifying the particle density as a Dirac delta function in order to reduce the domain of the integral 
to the particle’s worldline C. The orientation is chosen so that the spacetime volume is -kl = dt A dx A dy A dz. 


In (65), the curve C has constant x,y,z (due to the temporary choice of V) and / is a test function. Since the 
aim of this method is to induce the equation of motion of a particle from a fluid description, C is assumed to 
be an integral curve of V. 

In order to find the appropriate distributions for the particle versions of the magnetisation and polarisation, 
consider the following. Given a Killing 3-vector K € {i,j,k}, where i ■ i = 1, j ■ j = 1, k ■ k = 1 and 

i-j = i- k= j- k = 0,itis natural to introduce the de Rham current 


as 


(V-K)jfkl]=l^^l,-Kfdt, 

^For a summary of the conventions used in this Appendix, see Ref. m 


( 66 ) 
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where • represents the usual scalar product on 3-vectors and the subscript (v • Kj represents the distribution 
associated with the scalar V ■ K (likewise for other quantities). Expanding on the first of these equations: 

{V -k) j^il.-Kfdt 

= CD[fie-Kfdt] (67) 

and since Coictf] = {Cd A a) [/] for any 3-form a. Since the form degree of Cd is 3, Cd A dt = —dt A Co and 
hence 

1[/] = - (/7e • k) [dt A Cd) [/] ■ (68) 

Stripping off the test function and noting that in this case C = dt and ** 1 = —1: 

(iP • = - (/2e • i?) * (5 A Cd) ■ (69) 

Introducing the 4-vector /Xg = Hexdx + ^J■ey^y + ^J■ez^z, where ^ex is the a;-component of the vector /2e etc., it 
follows that the polarisation of the particle given as Vd = (v ■ Ka) dx^ is 

'Pd = {ic * Cd) fie, (70) 

for A S {1, 2,3} and similarly the magnetisation A4d of the particle is 

Md = {ic * Cd) fim- (71) 

From now on, the worldline C of the particle is taken to be arbitrary and /ig and fj,^ are 4-vectors orthogonal to 
the 4-velocity C of the particle. Furthermore, the timelike unit normalised 4-vector field F is C when evaluated 
over the image of C. 

The polarisation and magnetisation can be incorporated into a single 2-form, the polarisation 2-form 11 = 
C — F, where G is the excitation 2-form, given by: 

T\ = -V hV + ir{V ^M), (72) 

where V is the 4-vector describing the motion of the fluid. Using and 0, the distributional analogue of 
the polarisation 2-form can be written 


Ud — ★C'd a /ig — -k{'kCD A /Tin) 
= -kCD A /Tg — iy,^CD- 


(73) 

(74) 


J 


In order to establish a distributional analogue of the balance law (591, the distributional forms of the currents 
and must be formulated. The free current is found by noting that 


•free 


j 


= —qn * V, 


where q is the electron’s charge, and therefore via (65) we identify 

= qCd. 


(75) 


(76) 


Since the excitation form G satisfies d-kC = the field equation (62) shows that the bound current must 
satisfy j^ound _ Hence = —dkllD and the balance law (59) can be written 


dT^ = ^ [-qicix^Ff + * (df Mx^FA + * [Cx^F A *E) /} dr, 


(77) 


11 








where {X°‘} is a basis of translational Killing vectors and the 2-form E is defined as 


E — — Cl A /ie + *((11 A /Tin)- 

Further simplification occurs if we split df into its (fl-parallel and (ll-orthogonal pieces, leaving 


(78) 


• D 


{[ 


-qicix^F - E • {Cx^F) - Xc{{^x^F) ■ i^E)] / 




{a^df Mx-FA *e) I dr. (79) 

We now choose the stress distribution ansatz in order to satisfy this expression. Allowing T§ to be of the form 

Tg = -gin, X‘^)Cd + ic i^x^F A *E) A xCd. (80) 

where tt is a 4-momentum vector, with the second term designed to absorb the worldline-orthogonal pieces of 
(79), the stress balance law can be written in the simple form 

Vp [ix'^n F [ix'-F) ■ i^E) = -qi^ix-^F - E • [Cx-^F ), (81) 


where V is the Levi-Civita connection. Similarly, the spin balance equation (60), upon substitution of (80) 
becomes 


d<[/] = 


1 

2 1 


Cl“ (tt'’ + Hx^F) ■ i^E) + iix^F) ■ zx»E 

- o’" (7r“ -h {ix'^F) ■ ZpE) - [ix-^F) ■ ixbE fdr, 
and choosing the ansatz tr^ = for the spin de Rham current gives 

= -Cl“ (tt^ -h (ixbF) ■ i^E) - (ixoF) ■ ix^E 

+ (7r“ -h iix^^F) ■ i^E) + {ix<^F) ■ ix^E. 


(82) 


(83) 


A.l Substituting the Nakano-Tulczyjew momentum into the equations of motion 

Since the spin matching condition commonly used is the Nakano-Tulczyjew condition 0. it is logical to write 
the equations of motion (81) and (83) in terms of the 4-momentum P“ of the particle. This momentum may 
be expressed as 


P“(A) = -/' T^^Nb*N, (84) 

where T°‘^ is the stress-energy-momentum tensor and the 1-parameter family E\ of spacelike hypersurfaces is 
the set of leaves of a local foliation of spacetime with timelike unit normal N = Since the stress-energy- 
momentum tensor is related to the stress-energy-momentum forms T“ via T“ = *(T(A“,—)), note that 
for test 0-form /, 

[ r“AiV/= / *(r(A“,-)) AfV/. (85) 

Jm Jm 

Utilising the identity xa A N = xN A a where a is a 1-form and noting that the vector N is normalised as 
g{N, N) = —1, so that the volume form can be written *1 = —N A xN, this can be simplified to 



r“ A Nf 



iT{X‘^,N))N AxNf. 


( 86 ) 
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Since N = — , the integral can be split into a piece along the worldline C and another over the hyperplane 


Sa via 


r“ A TV/ = - 


IM 


IC 


dX 


* Nf 




IC 


\dX\- 


and stripping off the test forms yields the relation 

rSAN = -P^Cd a N, 

with TV“ = 1 ^. Using the stress forms given by (80) yields 

i-giiT, X‘^)Cd + ic (ix-FX *E) A *Cd) A P = -P'^Cd A P, 
which upon manipulation yields the relationship 


_ pa _ pab-p^ 


cb 


pc 

P-C] 


(87) 

( 88 ) 

(89) 

(90) 


between the momenta 7r“ and P“. Substitution of (901 into the equations of motion (81) and (83) yields the 
familiar equations 


d 

dr 

( P 


P“ + -r 

C^Pd 


1 , 


= -gP^^Cb - -E'’'=9“Pbc, 


dr 

found in the literature m- 


F^^E.dP'^ 

C-P, 


+ cM p“ + 


F’^’^E.dP'^ 

C-P, 


+ P'''=E “ - P“°E^ 


(91) 

(92) 


B Plasma wave electric field 

A sufficiently short and intense laser pulse propagating through a plasma drives a non-linear wave in the electron 
number density. For present purposes, the ions are essentially stationary over the timescales of interest because 
their charge to mass ratio is about three orders of magnitude lower than that of the electrons. 

Despite the recent focus on the three dimensional ‘bubble regime’ |26[ I27j , one-dimensional models remain 
useful for providing estimates, particularly in contexts such as the Stern-Gerlach force in the main body of the 
paper, where additional complexity leads to dramatically more difficult analysis. Some preliminary work on 
the subject of electrons in a one-dimensional maximum amplitude plasma wave has already been done in the 
context of non-linear electrodynamics |28j . and we use the same expressions for the plasma electron worldlines 
and overall electric field; a brief summary follows. 

Assuming the electric field is due only to the electron fluid and the ion background, the magnetic field 
vanishes, leaving only the electric field component in the direction of the propagation of the wave. The Faraday 
2 -form of such a wave driven by a laser pulse in the z-direction is thus simply 

F = E{^)dtAdz, (93) 

where ^ = z — vt is the phase of the wave. The electric field E is governed by the field and Lorentz force 
equations for a cold plasma: 

dE = 0, 

d 'k E — Qe^e ^ gion^ion ^ ^iont 
Vv = —ivF, 


m, 
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(94) 

(95) 

(96) 










where qion is the ion charge, Ue and riion are the proper number densities of the electrons and ions, 14 and l^on 
are the 4-vectors whose trajectories are the worldlines of the electrons and ions, respectively. For simplicity, the 
plasma electrons are assumed to be unpolarised and so the motion of the electrons is governed purely by the 
usual Lorentz force. 

For a plasma wave whose electron motion is much greater than the motion of the ions, it is reasonable to 
assume that the plasma ions are at rest, i.e. Ijon = dt with constant number density riion- A simple example is 
laser-plasma wakefield acceleration, where the plasma wave electrons oscillate far faster than the motion of the 
plasma ions. Seeking a 4-vector describing the motion of the electrons of the form 

K = HOdC - (97) 

the normalisation condition g{Ve,Ve) = —1 restricts ijjy giving 

Ve = i^iOdC - (98) 


where the negative sign is chosen so that the electrons move at speed less than the wave (for more details, see 

(99) 


]). In order to satisfy the Lorentz force (96), the electric field must be of the form 

, / 


E = 




QeT 


so that the function ly is similar to the electric potential. Similarly in order to satisfy the field equations (94) 
and (951 the electron number density can be written in terms of the ion density: 

^^7^9ion?^ion Zuion 


Qe 




T 


where Z is the degree of ionisation, Z = — , and iz must satisfy 

Qe 


d 


mz 


2ge7^ 


ly'^ — meZn\, 


{v^/. 


1/2 _ ^2 _ ^ 


= 0 . 


( 100 ) 


( 101 ) 


The square root in (101) imposes a lower bound on v. This results in a maximum amplitude oscillation for 
the plasma wave, and the maximum amplitude is known as the wave-breaking limit [29j . The wave-breaking 
limit may be obtained from (1011 by integrating from ^i, the minimum of v and hence a zero of E, to ^n, the 
maximum of E and turning point of v'. Since vi = !z(^i) = 7 and vn = 7^ (from (101)), it follows that 


F/max — ■\/2(7 1) 777.0.^U-ion- 


( 102 ) 


Given that the plasma wave attains its lowest possible value at ^i, integrating (101) between two subsequent 
zeros of the electric field and where 7z(^ni) = 7^(1 + gives 


to ; 


- meZuian (v\/- 7^ - ly + = 0, 


29274 

governing the electric field in the case of a maximum amplitude oscillation. 


(103) 
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